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previously determined effective-range parameters. Good description of the data is obtained with 
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alpha particles. 
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I. INTRODUCTION 



Nucleons in light nuclei have typical momenta that are small compared to the character- 
istic QCD scale of 1 GeV. At these low momenta, QCD can conveniently be represented by 
a hadronic theory containing all possible interactions consistent with the QCD symmetries. 
Effective Field Theory (EFT) provides a controlled framework for exploiting the separation 
of scales in nuclei. It is crucial to formulate a power counting that justifies a systematic 
truncation of the Lagrangian leading to observables with the desired accuracy. Nuclei offer 
a non-trivial challenge because one wants such a perturbative expansion in addition to the 
non-perturbative treatment of certain leading operators, which is required by the existence 
of shallow bound states. By now, mainly few-body systems have been studied within EFT, 
and, while much remains to be understood, many successes have been achieved [H, 0|. 

Similar to other approaches, the extension of EFTs to larger nuclei faces computational 
challenges 0,0]. As a first step in this extension, we specialized to very low energies where 
clusters of nucleons behave coherently 0, @, 0]. Even though many interesting issues of 
nuclear structure are by-passed, we can describe anomalously shallow ("halo" or "cluster") 
nuclei and some reactions of astrophysical interest. Since they are strongly bound, alpha 
particles play a central role in this framework. Many nuclear states have energies close 
to thresholds for break-up into alpha particles and nucleons, the most famous being the 
excited ("Hoyle") state of 12 C near the triple-alpha (3a) threshold. These states should 
be describable within the halo/cluster EFT, which is formulated with contact interactions 
among nucleon (N) and alpha-particle (a) fields. Together with the Na interaction, the 
aa interaction is an important input for such calculations. While we have studied the 
Na interaction elsewhere through both neutron-alpha (na) 0] and proton-alpha (pa) 0] 
scattering, we focus here on aa scattering. Consideration of this system is required before 
tackling other states with two or more alpha particles, such as 9 Be and 12 C. 

The internal alpha-particle dynamics is characterized by an intrinsic momentum scale 
Mhi associated with the binding mechanism. A naive guess is that this scale is set by the 
pion mass ~ 140 MeV. The aa interaction consists of the long-range photon exchange 
and short-range strong interactions. At low energies, the latter can be represented by con- 
tact interactions. The central issue is the relative importance of these contributions. The 
Coulomb interaction is non-perturbative for momenta smaller than around kc = ZiZ 2 a em iJ,, 
where a em = e 2 /4ii is the fine-structure constant, \i the reduced mass of the system, and 
Zj, i — 1, 2, the electromagnetic charge of the particles. Here fi = m a /2 and Z{ = Z a = 2 in 
terms of the mass and charge of the alpha particle, respectively, so kc ~ 60 MeV. At mo- 
menta much below 100 MeV, the deviation from pure-Coulomb aa scattering is dominated 
by the S wave [1, HJH3, 11, H, 13, 14, 15]. The large near-threshold ff-wave phase shift has 
been interpreted as resulting from a (J n , I) = (0 + , 0) state 14, 15, 16|, 17, 18] at an energy 
Er ~ 0.1 MeV above threshold in the center-of-mass frame, with a tiny width ~ 6 eV. 
The momentum corresponding to this 8 Be state sets a smaller scale M/ Q ~ Jcr = y / 2fiEn ~ 20 
MeV, which must arise from the larger underlying scale by a fine-tuning of the param- 
eters of the underlying theory. 

In the halo EFT, our goal is not to explain the mechanism of this fine-tuning, but instead 
to exploit its existence in order to describe a-cluster systems at low energies. We seek a 
description of these systems in an expansion in powers of the small ratio Mi a /Mhi- Power 
counting is dependent on how the various parameters scale with M/ G and Mm. The physics 
of the low-energy S state is conveniently discussed in the language of a dimeron field [l^] 



2 



with the quantum numbers of the low-energy composite state. This field is characterized in 
leading order by a fine-tuned mass A and a non-derivative coupling g to the aa state. In 
subleading orders more complicated couplings appear. 

It is not immediately obvious how the fine-tuned mass A relates to M/ Q . The simplest 
assumption is A ~ Mi Q M hi /n j20|. In the absence of Coulomb interactions, this reproduces 
the leading term in the effective-range expansion, and one has a shallow real or virtual 
bound state with a typical momentum ks ~ M/ Q . Strong interactions are non-perturbative 
for momenta of order ks and larger. Higher-order terms in the effective-range expansion 
appear as subleading corrections. This scenario is appropriate for S-wave NN scattering 



at momenta below Mh% ~ [2l|, |22j. For pp scattering, ks ~ 8 MeV and kc ~ 4 MeV. 



The Coulomb interaction can be included non-perturbatively in a straightforward way [23 
providing calculable contributions plus a renormalization of A/g 2 . 

The situation in aa scattering is somewhat different. The extremely low energy of the 
S'-wave resonance suggests that a smaller A might be necessary. An alternative fine-tuning 
assumes thus that A ~ Mf Q / y. Such scaling has already appeared in P-wave Na scattering 
jH|, 0] , and has striking consequences in S"-wave aa scattering. In the absence of the Coulomb 
interaction, the leading contribution for momenta k ~ M\ comes entirely from the unitarity 
term — ik in the inverse amplitude. To this order, the 8 Be system would be at the so-called 



unitary limit, exhibitin g co nformal invariance 2J, |25j, and 12 C would acquire an exact 
Efimov spectrum [I, [251. |26[| . As we will discuss later, this exact pattern is modified by the 
Coulomb interaction, which breaks scale invariance and thus moves the 8 Be ground state 
away from threshold. We will see that in leading order we can describe the ground state, 
and higher-energy, scattering data can be systematically accounted for in higher orders. At 
this point, we simply fit the EFT parameters to data. In the future, many-body methods 
jl, [i[ should allow the calculation of these parameters from the underlying EFTs, which in 
turn will have their parameters obtained from lattice QCD [27| . 

To the extent that the halo/cluster EFT is built on fields for clusters, it is closely related to 



phenomenological few-body cluster models [14|, [28j. The latter have had many applications, 
ranging from the general analysis of halo structures 29| to the successful confrontation 
with data for specific processes, such as decays of 12 C resonances 3(j. The emphasis here is 
instead on a systematic expansion of the most general interactions allowed by the symmetries 
of QCD. 

In this paper, we study the alternative expansion scenario and perform a detailed com- 
parison to phase-shift data for the aa system. As we will see, among the many findings our 
new method entails are a new expansion for the aa amplitude around the resonance, which 
differs from the effective-range expansion; the existence of inconsistencies between recent res- 
onance measurements and old scattering data; a more precise determination of low-energy 
parameters; and an extraordinary amount of fine-tuning. All these findings could have im- 
plication in the study of other alpha-cluster systems. The paper is organized as follows. In 
Sec. [Ill we briefly discuss the EFT formulation of the two-body system with Coulomb inter- 
actions, with some details relegated to App. [A] The proposed power counting is discussed 
in Sec. IIHI In Sec. IIV[ we present results of our fits to the + resonance parameters and aa 
phase shifts, and discuss the required fine-tunings. Our conclusions are presented in Sec. |V] 
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II. EFT WITH COULOMB INTERACTIONS 



We start with a summary of some basic ideas in halo EFT [1, 0, 0], extended to in- 
clude the main equations needed to deal with Coulomb interactions. These equations are 
a straightforward generalization of Kong and Ravndal's formalism (23| to include dimeron 
fields 1 . More details can be found in App. [Aj In Sect. II III we discuss the significant differ- 
ences in scales between alpha-cluster systems and the pp system considered by Kong and 
Ravndal 23 . 

We consider here the scattering of two alpha particles of charge Z a = 2 and reduced 
mass /i = m a /2, at a center-of-mass (CM) energy E = k 2 /2fi. We want to build an EFT 
that provides a controlled expansion for observables at momenta around the + resonance, 
k ~ k R = y/2nE R « 20 MeV ~ M lo < M hi . The energy of the resonance, E R ~ 0.1 
MeV, is much smaller than either the alpha-particle excitation energy E* a ~ 20 MeV or 
any energy set by pion exchange, such as two-pion exchange between alpha particles or 
one-pion exchange among nucleons, (2m 7T ) 2 /m a ~ m 2 ,/m^ m 20 MeV. Thus, we expect 
an alpha particle to behave in a first approximation as a rigid entity ("core" or "cluster"), 
which we represent by a scalar-isoscalar field <p. The smaller effects from deformations and 
other core-structure properties are accounted for in a derivative expansion. The breakdown 
momentum scale Mhi of this EFT will be set by the lowest-energy degrees of freedom that are 
not incorporated explicitly in the Lagrangian. Since these include the nucleons within the 
alpha particle — which can be resolved with a momentum ^m^E*, « 140 MeV — and the 
pions — which can be resolved with momenta of the order of the pion mass m n ~ 140 MeV- 
a reasonable estimate is Mh% ~ 



'm^E* a ~ m n ~ 140 MeV. The EFT provides an expansion 
of observables in powers of the ratio between the low-energy scales k and Mi and the high- 
energy scale Mhi- 2 Around the resonance k ~ k R , we expect an expansion parameter of the 
order of 1/7. As the energy increases the expansion deteriorates and one should not expect 
the EFT to be applicable to laboratory (LAB) energies much above Elab = k 2 / \i ~ 3 MeV, 
a conservative estimate that corresponds to a CM momentum of about 70 MeV. 

Since at low energies only S waves contribute significantly to aa scattering beyond pure 
Coulomb scattering, we introduce an auxiliary scalar dimeron field d with "residual mass" 
A, which couples to two a fields. (Observables are independent of the choice of fields.) We 
thus start with the following strong-interaction effective Lagrangian: 



C 



id + 



2m r 



a 



oft 



id + 



d 



fj 



dft<M> + (#) t d 



+ . . . 



where a = ±1 and g is a coupling constant. The ". . ." represent interactions with higher 
derivatives. The sign a is used to match the sign of the effective range tq. The existence of 



1 A similar extension was recenthz_considered by Ando and collaborators [32 

2 Recently it has been suggested [3l| that chiral symmetry is important for the properties of the 8 Be ground 
state. This argument is based on an extension of the halo EFT to include explicit pion fields. In fact, it 
had been pointed out before that the long-range interaction of alpha particles is related to two-pion 
exchange, but in the context of nucleoli scattering on an alpha particle built out of four nucleons. Since 
the pion mass is not much smaller than the typical nucleon momentum in the alpha particle, we see no 
rationale (for realistic values of the quark masses) for an expansion in powers of m^/Mhi in a theory with 
"elementary" alpha particles. 
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this sign merely reflects the non-physical, auxiliary character of the dimeron field |19l . |20 
In momentum space, associated with the displayed bilinear terms are the a propagator 



iS a (q ;q) = r— — — , (2) 

g - q 2 /2m a + le 



and the dimeron propagator 



iD d (q ;q) = — — — — — . (3) 

g - q 2 /^m a - A + %e 

Other terms in the Lagrangian can be treated as insertions of interaction vertices in 
Feynman diagrams or corrections to the propagators ([2]) and ([3]). 

Different assumptions about the dependence of A on M\ translate into different relative 
weights for the terms in the denominator of Eq. (j3D, where \q\ ~ y/2fiq = 0(Mi ). For 
example, in the present scenario all terms are of the same size, while in the NN case the 
energy and kinetic terms are of subleading order and the dimeron propagator reduces in 
leading order to iDd — > — ia/(A — ie). (The latter form gives rise in the amplitude not 
to a resonance but to a bound-state pole j20| at a momentum ks = i/ao, where ao is 
given by Eq. ({TBI .) In this section we illustrate the use of the Lagrangian (JTJ together 
with electromagnetic interactions to calculate aa scattering amplitudes. We postpone a 
discussion of the relative importance of various terms until the next section. 

Electromagnetic interactions are introduced into the effective Lagrangian ([!]) in the stan- 
dard way, that is, by both changing the derivatives into gauge-covariant ones, and intro- 
ducing gauge-invariant interactions generated by the electromagnetic field strength. For 
practical calculations one usually works in a fixed gauge. The Coulomb gauge is very con- 
venient, since it allows a clear separation between Coulomb and transverse photons. The 
former provides the leading electromagnetic interaction and is driven by the Sommerfeld 
parameter 

/ r \ Z (X em fl kc i . s 

vik) = — i — = Y ' ( 4 ) 

with a em = e 2 /47r the fine-structure constant and kc the inverse of the "Bohr radius" of 
the aa system. This parameter is enhanced by the presence of /x in the numerator, which, 
as we are going to see, makes t] large around the resonance and requires a non-perturbative 
resummation of Coulomb photons. On the other hand, transverse photons are suppressed 
flij by two powers of M\ j fi, and when in loops by extra powers of a em . Since numerically 
yU ~ 2Ml i /Mi , these interactions contribute to orders — where they can be accounted for 
in perturbation theory — beyond the precision we are working with. They can become 
significant only at energies comparable to the excitation energy of the a core. 

The scattering amplitude for two particles interacting via Coulomb plus a short-range 
interaction is given by j35[ (for more details, see App. |A|) : 

T = T c + T CSl (5) 

where Tq and Tcs are the pure-Coulomb and the Coulomb-distorted short-range amplitudes, 
respectively. Considering only S'-wave interactions, the Coulomb-distorted short-range am- 
plitude is parametrized in terms of the "Coulomb-corrected" phase shift Sq as 

2vr e 2iaa _ 2tt C 2 e 2ia ° 

CS ~ jfc(cot 5 C - i) ~ 1k c [K{rj) - H{rjj\ ' ^ 
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with the pure-Coulomb phase shift o"o given by Eq. (1A5I) in the Appendix and the Sommerfeld 
factor 

,2 27r V 



CI 



The H function is given by Eq. (1A15I) . For real 77 it can be expressed as 

H{rf) = Re ty>(l + ir,)} - In 77 + 

2,7] 



(7) 



(8) 



in terms of the digamma function ip{z) = (d/dz) In T(z). The other (real) term in the 
denominator of Eq. (jHj) is the Landau-Smorodinsky K function [371 ] . 



JiC(i7) = ^ (cot + 



(9) 



At low energies it reduces to the effective-range expansion (ERE) in the presence of Coulomb 
interactions, 



K( v ) 



2k 



c 



a 2 



4 



(10) 



where a , r , Vo, ■ ■ ., are the scattering length, effective range, shape, . . ., parameters. 
As shown in App. [Aj the calculation using the EFT Lagrangian ([!]) leads to 



T C s = --tfe 2 ^ 



a s (J^r^ k - 2k c H{r]) 



/j, 2 g 2 



(11) 



where is the Coulomb-renormalized mass parameter of the EFT Lagrangian, 



= A(k) 



a- 



2vr 1 D - 3 



2k 



o 



D — A 



In 



2k 



- 1 



c 



3 -c E 

2 



(12) 



with k the renormalization scale and D the dimension of spacetime. Equation (TTTj) is in 
the form of the Coulomb-modified ERE with scattering length and effective range given 
respectively by 

(13) 



do = -v- 



r = -cr 



2tt 
H 2 g 2 



2ttA( r ) ' 

From Eq. ffTTl) it is clear that the effect of a non-perturbative Coulomb dressing of the 
strong-interaction amplitude — apart from multiplying the amplitude by C^e 2tcr ° and from 
renormalizing the short-range parameters — is to effectively replace the unitarity term —ik 
by —2kcH{rj). In order to estimate the relative sizes of the various contributions to Eq. 
f TTTj) . we turn now to a discussion of power counting. As we are going to see, at each order 
the EFT in the two-body system is equivalent to a truncation of the Coulomb-modified 
ERE. 



III. POWER COUNTING 

In this section we elaborate on the proposed power-counting scenario mentioned in the 
Introduction. We also give the aa amplitudes that are used in the next section for a 
comparison to phase-shift data. 
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In the present strongly-tuned scenario one considers A ~ Mf / 'n and g 2 /In ~ M^/// 2 , 
with other couplings scaling with Mhi according to naive dimensional analysis. In this case 
the contribution of the bare dimeron propagator to the scattering amplitude at momentum 
k ~ Mi comes not only from the dimeron mass, but also from its kinetic term. The sim- 
plest strong-interaction contribution is a bare dimeron propagator attached to initial and 
final external legs, which (as seen by multiplying Eq. by g 2 ) is of 0{2nM hi / fiM 2 ). A 
non-Coulomb bubble diagram times another dimeron propagator brings an extra factor of 
0(fig 2 k/2nA) = 0(Mhi/Mi ). As a consequence, the bubble-chain resummation is neces- 
sary. The resulting denominator acquires the form (— 1/ao + r^k 2 /2 — ik). Interesting in 
this power counting is the fact that at momenta k ~ M\ the first two terms are suppressed 
by Mi /Mhi compared to the last one. Therefore, all that is left at LO if the Coulomb 
interaction is turned off is the unitarity term l/(—ik). In this limit, the 8 Be system exhibits 



conformal invariance [24, 251 and the corresponding three-body system, 12 C, acquires an 
exact Efimov spectrum 2j, |25|, [26] . This scenario is a possible realization of the unitary 
limit. 

This picture is significantly modified when the Coulomb interaction is turned on. In a 
non-relativistic system an 1/r potential breaks scale invariance and introduces the scale kc in 
the propagation of two charged particles. As we have seen in Eq. (TTTT) . the unitarity term is 
modified. The balance between strong-interaction terms and Coulomb-modified propagation 
depends on both the intrinsic strong-interaction scale and kc- We note that the scales are 
very different in aa than in pp scattering. While for the latter kc/ks ~ 1/2, for the aa 
system, kc = a em Z 2 m a /2 ~ 60 MeV and kn = ^m a Efj ~ 20 MeV, so kc/kn ~ 3. For 
momenta k ~ kR, we are therefore in the deep non-perturbative Coulomb region. 3 This case 
corresponds to large values of rj, where the function 2kcH{rj) is significantly different from 
the usual unitarity term ik. Using Stirling's series [36| . 



lnr(l + 2 ) = il ll2 , + ( 2 + i)h 1 ,-, + I i-- s L 5 + ..., (14) 
in Eq. (jSJ) gives 

HM -j^ + mf + --- + 7^- 1 - (15) 

The unitarity term is thus replaced by 2kcH{rj) ~ k 2 /Qkc- This term is now a factor k/Qkc 
smaller in magnitude than the unitarity term in the absence of Coulomb, and comparable 
to the effective-range term coming from the dimeron kinetic term. This can be captured 
automatically if we take 3kc ~ M^, as it appears to be the case numerically. 

If all higher-order parameters are natural, in the S wave the shape-parameter term 
Vok 4 = 0{MfjM^ is down by (Mi /M hi ) 2 compared to the effective-range term r k 2 = 
0(M 2 /Mhi), and so should be each of the successive terms. On the other hand, the D- 
wave scattering- "length" term a 2 k 4 = 0{Mf /Mf ii ) is down by (Mi /Mhi) 6 compared to 
the S- wave scattering length ao = 0(Mhi/M 2 ), and further energy dependence and higher 
angular momenta bring in further powers of (Mi /Mhi) 2 . 

In this power counting, the Coulomb-distorted short-range amplitude T C s is given in LO 



This fact suggests that one might be able to develop a perturbation scheme in powers of k^/kc 
y. Z 2 



2^/E R /m a /a em Zl 
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by Eq. ([TT]) . Including corrections in perturbation theory, we obtain up to NLO: 



rp _ 2-7T 2 2j<r 

J as — % e 

/i 



■l/a + k 2 r /2-2k c H 4 (-l/a + P r /2-2k c Hy 



(16) 



where "Po is given by a higher-derivative term in the Lagrangian (CQ). This corresponds to 
an expansion of the ERE formula, 

T (ERE) = _fn_ U^e 

cs fi -l/a + k 2 r /2-k 4 Vo/4-2k c H ' 1 ; 

for small P - To this order, all waves higher than S are purely Coulombic. Higher orders 
can be calculated similarly 

Equation (Tl6l) holds for generic momenta k ~ M/ Q . However, it fails in the immediate 
proximity of k^. This situation is familiar from another application of the halo EFT to a 
resonance The power counting works for k ~ Mi Q except in the narrow region \k — k^\ = 
0(Mf / Mhi) where the LO denominator nearly vanishes and a resummation of the NLO term, 
here associated with the shape parameter, is required. As one gets closer to the resonance 
momentum kf>, higher-order terms in the ERE are kinematically fine-tuned as well. This 
happens because the imaginary part of the denominator is exponentially suppressed by a 
factor exp(— 27r?7^) ~ and the real part gets arbitrarily small. (For na scattering, 

multiple kinematical fine-tunings are prevented in the P3/2 wave by the presence of an 
imaginary term ik 3 ~ Mf Q [H, 0].) 

This kinematical fine-tuning is not a conceptual problem. From the EFT point of view, 
each new fine-tuning is equivalent to reshuffling the series and redefining the pole position. 
Such a procedure works fine with a small number of kinematical fine-tunings, but is not 
practical in the present case. A better alternative is to perform an expansion around the 
resonance pole position. The situation here is analogous to the NN system, where we can 
choose to expand the amplitude around the bound-state pole [22| rather than around zero 
energy. 

A great simplification results from the fact that the resonance lies in the deep Coulomb 
regime, where Eq. (j!5p provides an accurate representation of H up to the precision we 
are considering. The real terms shown in Eq. (fT51) can be seen as an expansion in powers 
of ~ (k/3kc) 2 = 0{k 2 /M^). Of course in an asymptotic expansion at some point the 
remaining terms can no longer be expanded; at that point the remainder should be treated 
exactly. In lowest orders, however, we can use the successive terms shown in Eq. (TL51) : 
numerically, the terms up to work to better than 3% for Elab = 3 MeV. 

The expansion ([151) makes the physics around the resonance quite transparent. Since 
the "size" of the resonance, 1/kpt, is much larger than the Bohr radius 1/kc, the Coulomb 
interaction is effectively short ranged, and the real part of H is ERE-like. This expansion 
matches well with the expansion for the Landau-Smorodinski function, which is in pow- 
ers of (k/Mhi) 2 . In Tcs, n ot only the k 2 terms, but also higher-order terms from strong 
and Coulomb interactions have comparable sizes. We can thus lump these terms together, 
defining 

n, = r„-J-, P = n + J^, (18) 

and so on. 
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In this case we can rewrite up to NLO: 



T, 



2tt 



cs 



nl 2ia 



H -1/ao + r k 2 /2 - Vok^/i - ikC 2 

2tt gg° 

7 r (P - k 2 R )/2-V (k 4 -k 4 R )/4-tkC 2 



2tt 



£(2 e 2io- 



(A' 



4 -4) 



f (A; 2 - k 2 R )/2-ikC 2 4 {r {k 2 -k 2 R )/2-ikQ 



2)2 



+ 



(19) 



LO term 



where 



a r 



Vo 
a fo 



NLO correction 



+ 



(20) 



LO term NLO correction 

From this expression one sees directly that, indeed, k R ~ Mi , with corrections of 
0(Mf / M^j)- Note that we keep the exact form of the imaginary term in Eq. (TT9]) : even 
though it is negligible at k ~ k R , it has an important exponential dependence on the energy 
responsible for keeping the phase shifts real in the elastic regime. 

When a < and f < 0, and thus r < l/3k c , we have k\ > and the two poles of 
Eq. (jl~9l) are located in the lower half of the complex-momentum plane very near the real 
befits a very narrow resonance. The amplitude Tcs can be written in terms of the 
resonance energy Er = k R /2fi and the resonance width T(E) as 



Tcs 



2ire 2ia ° 



T{E)/2 



fiy/2fIE E- E R + iY{E)/2 



(21) 



One finds that 



T(E) = T(E R ) 



e 2irk c /k R _ -y 



2nk 



c 



kc/ku 



■(E-Ef 



(22) 



LO term 



NLO correction 



where 



Auk 



c 



+ 



^0 



+ 



LO term NLO correction 



(23) 



The width is very small because of the large value of 2nkc/kR in the exponential. 

In the form of Eqs. (f2~T|) and (122]) we can keep Er and F(Er) fixed at each order in the 
expansion. Note that these equations do not change to this order if one makes a different 
choice — e.g., {k 2 — k 2 R ) 2 instead of fc 4 — k R — for the form of the Vq term in Eq. (jT9l) . Since 
at the resonance, Sq(Er) = ir/2 and [(16q(E) /cIE]e r = 2/T(Er), the behavior of the phase 
shift around Er is fixed. We turn now to a test of the expansion (1T91) mandated by our 
power counting. 
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IV. COMPARISON TO DATA 



In this section we briefly describe the experimental situation regarding aa scattering at 
low energies. Our predictions are then compared against the available data and discussed 
in detail, with a particular emphasis on the scalings of the EFT parameters. 

Scattering data at low energies are not abundant. At energies below E L ab = 3 MeV 
data were obtained, and a 
that we were able to find 



)hase-shift analysis performed, in Ref. [8j. All later references 
LH . 15 1 start at higher energies. For example, Ref. 0] covers 
the region E LAB = 3-6 MeV, Ref. E LAB = 3.8-12 MeV, and Ref. [U E LAB = 5- 
9 MeV. These data show that aa scattering at Elab ~ 6 MeV is dominated by the S 
wave, thanks to the presence of a (J\ I) = (0 + , 0) resonance immediately above threshold, 
identified as the 8 Be ground state 15[ . Early determinations of the + resonance energy were 
performed in reactions like n B + p — > 2a + a [16|. Later measurements and careful analysis 
of the scattering of 4 He atoms off 4 He + ions |l7l. 18| improved the determination of the 
resonance energy and width to their currently accepted values, E^ AB = 184.15±0.07 keV and 

^lab = 11-14 ±0.50 eV. The resonance CM momentum is thus k B = y fiE^ AB « 18.5 MeV. 
At Elab ~ 6 MeV the .D-wave phase shift crosses tt/2, indicating the position of a broader 
resonance associated with the first excited, (2 + , 0) state [l5j |. The G wave does not become 
comparable to the D wave until Elab ~ 20 MeV, interpreted as the region of an even 



broader (4 + , 0) state [15 



resonance 



energy 
10 3 fm, 



The low-energy data can be studied with the ERE. When used together with the 
16| . the data from Refs. jl, 0] give for these parameters 

1.096 fm, and Vq = —1.654 fm 3 . They provide a good descrip- 
6 MeV, which lends some credence to these 



do = —1.76 x lir tm, Tq = i.uyu mi, anu ro 
tion of the available phase shifts up to Elab 
numbers. However, it was pointed out 12| that without input from the + resonance width 
these parameters have large uncertainties. Inclusion of both resonance energy and width 



from Ref. 17J] reduces these uncertainties considerably [13J: ao = (—1.65 ± 0.17) x 10 3 fm 



r = (1.084 ± 0.011) fm, and V 
from Ref. 



1.76 ± 0.22) fm 3 . Since the later, more precise data 
18j are consistent with Ref. [13], the ERE parameters from Ref. [l3| can be seen 



as a reasonable parametrization of the existing data. 



Since it is well-determined 



Here we use the phase shifts compiled in Table II of Ref. [141 
experimentally, and due to its relevance to the triple-alpha process, we use the + state as an 
important constraint. This is in line with the EFT approach, where lower-energy observables 
have smaller theoretical errors. It provides a relationship among our EFT parameters and, 
consequently, reduces the number of variables to be adjusted at each order in the power 
counting. Below, we also use the ERE from Ref. [l3[ for orientation, and comment on the 
extremely large value of the scattering length do, which suggests a large amount of fine- 
tuning in the parameters of the underlying theory away from the naturalness assumption. 
In contrast, both r^ 1 ~ 180 MeV and Vq 1 ^ 3 ~ 170 MeV correspond to natural scales 
comparable to the pion mass. At this level there is no evidence against our initial estimate 
of the expansion parameter Mi /Mhi ~ 1/7. Note that here we do not include dimeron fields 
for resonances beyond the ground state, and therefore cannot go beyond the energy region 
where the D-w&ve resonance is significant. 

In the power counting we are proposing for the aa system, the amplitude Tcs for generic 
momenta is given up to NLO by Eq. ([19]) . As previously discussed, this expression combines 
the deep-non-perturbative Coulomb approximation (j!5p for the function H with the expan- 
sion around the resonance pole, which avoids the need for multiple kinematical fine-tunings. 
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FIG. 1: Left panel: S- wave phase shift <5g as a function of the laboratory energy Elab- The EFT 
results in LO and NLO are given by the (blue) dotted and (red) solid lines, respectively. Our ERE 
fit using Eq. (fTT|) is given by the dash-dotted line and the empirical phase shifts [14|] as (black) solid 
circles with error bars. Right panel: analogous results for the Landau-Smorodinsky ^-function. 



In LO, the two parameters cto and f can be obtained from a fit to the resonance position 
and width. At NLO, scattering data are needed to determine Vq. 

Figure [H shows the results of our fit to the available S'-wave phase shifts below Elab = 
3 MeV, including the resonance position and width. The latter control the steep rise of the 
phase shift at very low energies. In the region above the resonance, where scattering data 
are shown, the LO curve is a prediction, which is consistent with the first few points but 
then moves away from the data. The NLO curve has an extra parameter, which here was 
determined from a global x 2_n t to scattering data shown. As expected from a convergent 
expansion, the description of the low-energy data improves with increasing order. At about 
3 MeV and above, higher-order contributions are expected to be significant, as suggested 
by the discussion on the relevant scales and manifest in the growing difference between the 
NLO curve and both LO curve and data points. Also shown are results from a fit using the 
conventional ERE formula, Eq. (j!7p . in order to stress the differences between this and our 
EFT approach. The ERE formula includes some of the contributions of higher order in the 
EFT. 

In Table [J we give the values of ERE parameters that we extract from the fits in Fig. [I], 



and compare them with the values [13| obtained from effective range theory. At LO, ao and 



ro come out consistent with the values given in Ref . [13J , r$ in particular. The changes due to 
NLO corrections, however, worsen this initial LO agreement. The reason for this deviation 
is most likely due to the calculation of the width constraint in Ref. [13]. Its Eq. (4) reads 

dh 1 7i 1 / 2 \ /^>,\ 

dk^ R) - « e^o/ kR - 1 = W c ( r ° " V ° kR ) ' (24) 

where h(r]) = ~Re[H(r])]. Following this reference's procedure, we were able to reproduce the 
quoted value of the width (6.4 eV) only when dh(r]ji)/dk 2 was approximated by 1/12/^ [see 
Eq. (fT5l) ]. This amounts to ignoring the electromagnetic piece of Pq in Eq. (|23|) . which is 
inconsistent since the strong piece contributes at the same order. It explains why the value 



for ro in Ref. [13J agrees quite well with ours at LO but disagrees at NLO. With this larger 
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fo, one can also understand why Ref. [13j obtains a smaller ao (see the following discussion 
on the scaling of ao). When we repeated Ref. Hf's procedure including the width constraint 
consistently, we obtained essentially the same values as in our EFT fits. This updated ERE 
fit is also shown in Table [B 

In agreement with our expectations, the effective range shows a natural size, r ~ 
1/(180 MeV). The shape parameter also has a natural size Vo ~ 1/(170 MeV) 3 at NLO 
which is in agreement with our a priori estimate. The relative errors in ao and f at LO 
are estimated to be of the order of the EFT parameter expansion, M^/M^ ~ 1/7 ~ 15%. 
At NLO, they are essentially given by the precision of the most recent measurement of the 
resonance width 18j], which lies between 4-5%. The uncertainty in Vo, given by the x 2_ fiL 
is of the same order. Note that the small relative error in ro compared to the one in fo is 
due to the former being an order of magnitude larger than the latter, as we discuss below. 
The NLO errors found here are a factor of two smaller than the ones obtained by Ref. 13 . 

One should stress that the accurate value of the resonance width, T R = 5.57 ± 0.25 eV, 
imposes tight constraints on our fits, through a and r . There is a significant improvement 
in our NLO fit and overall agreement with data, but the theoretical curve is not able to 
cross the error bands of all scattering points below 3 MeV, as reflected in a x 2 /datum ~ 4. 
In principle, a better agreement should be achieved by an N 2 LO calculation. However, that 
introduces an extra parameter that is mostly determined by the scattering data, and this 
agreement could mask any possible inconsistencies that the phase shifts might have with the 
resonance parameters. The high NLO x 2 /datum suggests that the resonance width and the 
5*- wave scattering data set are not compatible with each other or, at least, one of them has 
overestimated precision. As a test we performed fits to scattering data without the input 
from the resonance width, within our NLO EFT and the conventional ERE framework. In 
both cases, description of S-wave phase shifts is much better but the resonance width is 
underpredicted, T R = 4.9 ± 0.6 eV with ERE and T R = 2.87 ± 0.23 eV with EFT. The 
ERE result is still consistent with the measured T R thanks to its large error bar. In EFT, 
where lower-energy data have higher priority, the discrepancy is amplified. The problem is 
even more pronounced if the fit is performed using data up to 2.5 MeV instead of 3 MeV: 
the results T R = 4.2 ± 0.6 eV with ERE and T R = 2.93 ± 0.34 eV with EFT fall beyond 
the quoted experimental error bars. Oddly, this tendency continues as one lowers the upper 
limit in the fit. Reanalyses of the existing low-energy data or even new measurements seem 
necessary to resolve this discrepancy. 

An astonishing feature of the aa system is the very large magnitude of ao, even if com- 
pared to the large scattering length observed in the two-nucleon system. The latter is 
evidence of a fine-tuning in the QCD parameters, which gives rise to an anomalously low 





a (10 3 fm) 


r (fm) 


Vo (fm 3 ) 


LO 


-1.80 


1.083 




NLO 


-1.92 ±0.09 


1.098 ± 0.005 


-1.46 ±0.08 


ERE (our fit) 


-1.92 ±0.09 


1.099 ± 0.005 


-1.62 ±0.08 


ERE [13] 


-1.65 ±0.17 


1.084 ±0.011 


-1.76 ±0.22 



TABLE I: ERE parameters extracted from EFT fits in the first two orders, compared with values 
from two ERE fits, our own and Ref. [13j's. 
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momentum scale. (This fine-tuning can be seen as the proximity of the observed pion mass 



to the critical value ~ 200 MeV where the NN bound states have zero energy [27|, |38| and 



the triton displays the Efimov effect .39].) Certainly this fine-tuning scale propagates to 



heavier systems. However, the enormous value of a in the aa system is suggestive of a 
more dramatic fine-tuning, with electromagnetic interactions playing a crucial role. 

The parameters of the + resonance are indeed quite surprising. As we remarked, this 
resonance is associated with two poles of Tqs, the momenta of which are much smaller in 
magnitude than kc- As a consequence, there is an exponential suppression of the width, 
which is evident in Eq. (123ft . Yet, the width is somewhat large given the position of the 
resonance. It is remarkable that in order to fit both the resonance position and width one 
needs fo — —0.13 fm, nearly a factor 10 smaller in magnitude than the Coulomb contribution 
to the energy-dependent term, — l/3kc — — 1.2 fm. That is, although both r ~ 1 fm and 
1/3 fee ~ 1 fm as expected from dimensional analysis, they approximately cancel producing 
a balance that is about 10% of either. This balance would itself have natural size if the 
width were a factor of 10 smaller than measured, or else if Jcr were 15% larger. Instead, the 
existing cancellation effectively leads to either |f | ~ or |f | ~ Mi /M^. 

Since the position of the resonance determines to a first approximation aofo, the seemingly 
accidental cancellation in fo translates into an exceptionally large do, which is effectively 
|ao| ~ fi>/Mfo or |ao| ~ Mfa/Mfo. The Coulomb-modified parameter is thus suppressed 
by another power of Mi Q /Mhi compared to the assumed scaling when Coulomb interactions 
are turned off. 

The values of ERE parameters show that, as anticipated by the power counting we 
proposed, the strong-interaction effective-range term is comparable to both the scattering- 
length term and the electromagnetic H . A power counting based on the scalings appropriate 
for AW scattering would instead miss the resonance in LO unless one also counted 3kc as 
Mio, which is a numerical stretch. Even so, the width would be far off, since l/3kc 3> |fo|, 
and convergence questionable. 

Although the results presented in this section were based on Eq. f|T9|) . and thus assumed 
the validity of the expansion (fT5l) . we have checked that no qualitative changes arise if one 
bases the EFT fits on Eq. ffTBl) instead. In the latter case, one is simply retaining at each 
order some terms that truly belong to higher orders, at the cost of much less transparent 
analytical expressions. 



V. CONCLUSION 

In this paper, we have studied aa interactions within the framework of halo/cluster EFT. 
From this perspective, the aa system is most extraordinary. 

EFT is model independent to the extent that includes all interactions allowed by known 
symmetries. It is equivalent in the aa system to a truncation of the ERE. The halo EFT 
can be viewed as a formulation of the ERE that can be used for other systems made of 
shallow-bound alpha particles and, together with results from Refs. p, H, 0], also nucleons. 

We have designed a power-counting scheme that seems realistic for the aa system, but 
requires even more fine-tuning than the one in the NN system. It assumes that the param- 
eter A in the EFT Lagrangian scales as Mf / fi in the absence of the Coulomb interaction. 
Interesting in this scenario is the fact that at LO without Coulomb the 8 Be system would 
exhibit conformal invariance and 12 C would display an exact Efimov spectrum. These exact 
features are broken by the Coulomb interaction but some remnants of this behavior are 
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manifest in the experimental spectra, such as the shallowness of the + 8 Be resonance. 

We have incorporated two new ideas in the halo/cluster EFT. The most effective way of 
implementing an EFT with resonances seems to be an expansion around the resonance pole, 
since it avoids the need for multiple kinematical fine-tunings and improves convergence. Our 
study was significantly simplified by treating the Coulomb function H(rj) in a power series 
expansion in k/kc- This simplification may be very relevant when treating interactions with 
several alpha-particle clusters. 

We made use of the precise measurements of the + resonance properties [l8|, together 
with the (rather old) existing S-w&ve scattering phase shifts [s| in order to produce LO and 



NLO fits that fix our ERE parameters. The result is in slight disagreement with Ref. [13 
probably due to an approximation made in its calculation of the width. The uncertainties 
in our ERE parameters are smaller by a factor of two in comparison with this reference. A 
systematic improvement is seen in the theoretical phase shifts, and reasonably good overall 
description is obtained at NLO. Yet, the rather high \ 2 when the resonance width is fitted 
and the deterioration of the calculated width as one restricts the fit to lower-energy scattering 
data suggest that either the resonance width or the S*-wave phase shifts have overestimated 
precision. 

Despite the phenomenological success of the proposed power counting, the necessary 
cancellations between strong and Coulomb interactions, which seem accidental, are puzzling. 
In the absence of Coulomb interactions, the regulator-dependent part of the bare dimeron 
mass A is 0(M|j//i). As we can see from Eq. (|T2l) . the electromagnetic ^-dependent pieces 
do not change this expectation, since they are ~ 2kcfig 2 /2n = 0(M^ i / jj). Therefore, barring 
cancellations between the finite pieces of A and loops, we would expect the renormalized 
mass A^) to have the same size, which would set a scale for |a | at 0(1/ Mm) ~ 1/2&C ~ 2 
fm. Yet, the observed resonance energy (together with a natural-sized effective range rn) 
tells us that A^ = 0(Mf / jx). A similar scaling for A^ emerges in na scattering [5|, |6j, 
but here we need a further cancellation between the finite pieces in A and Coulomb loops. 
This brings |ao| up by two orders of magnitude, or more than two inverse powers of the 
expansion parameter M io /M hi ~ 1/7, to 0(M hi /M? ) ~ 100 fm. 

Uncomfortable as this fine-tuning by a factor of nearly 100 in the energy-independent 
part of the amplitude Tcs might be, it is not the whole story. As we have shown, both 
the resonance width and the higher-energy phase shifts require a fine-tuning also in the 
energy dependence of the amplitude. This ~ 90% additional cancellation (together with the 
observed resonance energy) further enhances |<xo| by another factor of 10, or about another 
inverse power of Mi /M h i, leading effectively to 0(MlJM? ) ~ 700 fm, which is indeed 
the order of magnitude we obtain in our fit. This fine-tuning of a factor of ~ 1000 in aa 
completely overshadows the fine-tuning of ~ 10 — from (or from the pion decay constant 
f n ~ 92 MeV) down to l/|a | ~ 8 MeV — in the NN 1 S channel. It has important 
consequences: for example, if the strong-interaction effective range r were just 15-20% 
larger it would make the ground state of 8 Be bound, presumably with far-reaching effects in 
nucleosynthesis. 

This context frames a fascinating picture for the aa system: a nearly conformally invari- 
ant system that is plagued by cancellations between strong and Coulomb interactions. Of 
course, fine-tuning has long been discussed in nuclear physics, but usually in connection to 
the position of the Hoyle state of 12 C (see, for example, Ref. jioj). An immediate question 
is to which extent the latter arises from the fine-tunings we discussed here. 

Together with previous and ongoing work on the na [1, @] and pa systems 0] , our results 
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provide a framework for the description of cluster states in nuclei using EFT, which can 
be seen as a generalization of the ERE to systems that involve more than two bodies. An 
important extension of this work is to systems of more than two alpha particles. The strong 
aa interaction obtained here, together with an exact treatment of the Coulomb interaction, 
should provide the necessary ingredients. Since k c /kR ~ 3 as discussed above, one might 
be able to simplify the calculation of Coulomb effects by developing a perturbation scheme 
in powers of k^/kc from the start. Moreover, we conjecture that the 12 C Hoyle state is 
a remnant of an Effmov state that appears in the unitary limit. More complex a-cluster 
systems could also be studied with the EFT presented here. 
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APPENDIX A: COULOMB GREEN'S FUNCTION AND EFT 



The scattering amplitude for two particles in their center-of-mass (CM) system interacting 
via Coulomb and a short-range interaction is given by (35| 



T = T, 



c 



T cs = (xi'^Vclk) + ( x [J>\Vs\^n 



(-)i 



■,(+)\ 



(Al) 



where \k), an d ) represent free and (incoming/outgoing) states of momen- 

tum k for pure-Coulomb and Coulomb-distorted short-range interactions, respectively, while 
Vc (Vs) is the Coulomb (short-range) interaction operator. In coordinate space, the Coulomb 
wave functions can be written as 23 



(rlxjb ) = xC } ( r ) = e ' 2 r(l ± irj) M{=fir], 1; ±ikr — ik ■ r) e 



ikr 



(A2) 



where M(a, b; z) is the Kummer function. From M(a, b; 0) = 1 [36|, one obtains the impor 
tant properties 

2irri 



*(0) X? } (0) = eT(l =F <i7)r(l ± i V ) 



X? ) *(0)xl ±) (0) = e--T(l±^ 



q2 e ±2iao 



(A3) 
(A4) 



where <To is the Coulomb phase shift for the partial wave 1 = 0. The general expression for 
the Coulomb phase shift, 



oi = argT(/ + l + irf) 



2i 



In 



r(Z + 1 + irj) 
V{1 + 1 - irj) 



(A5) 
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is defined from the partial-wave expansion of the pure-Coulomb amplitude 

,2i* t _ 1) 2 7T 



The explicit solution 



p 1=0 



2ik 



P/fCOSI 



fc{9) . 



(A6) 



fc{0) 



2k 



o 



csc 2 #/2exp irj ln(csc 2 9/2) + 2iao 



(A7) 
I 2 , which 



leads to the well-known Mott scattering cross section, au = \fc(@) + fc(n ~ 
holds at very low energies. 

The main ingredient in the calculation of the Coulomb-distorted short-range amplitude 
is the Coulomb Green's function at energy E, 



G$\E) 



1 



E- H c ±ie 



2/i 



d 3 Q \X^)(X^\ 
(2vr) 3 2/iE - q 2 ± ie 



(A8) 



Using the Lippmann-Schwinger equation, one is able to express l^j^) in terms of multiple 
insertions of the operator Gq^Vs acting on the Coulomb states |xi ± ' > ) 23], allowing the 
amplitude to be written as the sum 



Tcs = Uxl- ) \Vs(G¥v s ) n \xl +) ) 



(A9) 



n=0 



This amplitude can be decomposed into partial waves, 



T C s = ^TiPi(cost 



(A10) 



in terms of the "Coulomb-corrected" phase shifts 8f. 

If the deviation from the Mott cross section is mostly 5" wave, then the calculation of Tcs 
proceeds from the Lagrangian ([1]) similarly to Ref. [23j. The purely strong matrix element 
(p'\Vs\p) is defined as the amplitude td = g 2 Dd(E; 0) for aa scattering via a bare dimeron 
two-point function, with the external legs amputated, evaluated in the CM: 



(p'\V s \p)=g 2 D d (E;0) 



ag< 



E - A + ie 

Explicit evaluation of Eq. (TAH using Eqs. (lA"3l) . (IA"4|) . and (lAlip then gives 

\ 2 



(All) 



T, 



cs 



° 9 r 2 2ia 

E - A + ie v 



1 



E- A + ie 



- ME) + 



E- A + ie 



-ME) 
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This resummation of the Coulomb exchanges is illustrated in Fig. [21 It leads to 

-1 



T< 



cs 
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FIG. 2: Resummation for the LO expression of Tcs, represented by the (red) shaded rectangle. The 
double line stands for the dimeron propagator. The (blue) shaded ellipse represents the propagation 
in the presence of Coulomb photons (dashed line) but absence of a dimeron. 



where Jq is given by 



/i I K 

2tt I D - 3 



+ 2k 
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(A14) 



with D the dimensionality of spacetime, k the renormalization scale, kc = Z 2 a em fi = kr), 
Ce = 0.577... the Euler-Mascheroni constant, and 



H{rf) = ip{irf) + hi{irf) 

Jut 1 1 



in terms of the digamma function ip, which obeys [3J 

Re [ip(ir])] = Re [0(1 + ir/)] , Im [i/j(ir])] 



1 7T 

2^ 2 
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The dimeron mass gets renormalized by the non-perturbative free-particle and Coulomb 
loops, 



27rA^) 2vrA( 
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and Tqs finally becomes 
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